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ABSTRACT 

There exists an upper limit on the mass of black holes when the cosmological 
constant A is positive. We study the collision of two black holes whose total mass 
exceeds this limit. Our investigation is based on a recently discovered exact solution 
describing the collision of Q = M black holes with A > 0. The global structure of 
this solution is analyzed. We find that if the total mass is less than the extremal 
limit, then the black holes coalesce. If it is greater, then a naked singularity forms 
to the future of a Cauchy horizon. However, the horizon is not smooth. Generically, 
there is a mild curvature singularity, which still allows geodesies to be extended. 
The implications of these results for cosmic censorship are discussed. 



1. Introduction 



It is almost 25 years since Penrose first proposed his cosmic censorship con- 
jecture [1]. Roughly speaking, this conjecture states that naked singularities do 
not form from realistic physical processes. This conjecture is widely believed to 
be true and has become the cornerstone of our understanding of gravitational col- 
lapse and black holes. But despite extensive work over the years, we are still far 
from having a proof. One appears to need global existence results for strong-field, 
nonsymmetric solutions, which are extremely difficult to obtain. 

When a proof seems difficult, it may be easier to obtain a convincing counterex- 
ample and show that the conjecture is false. Many attempts have been made. One 
of the earliest is the following. Static charged black holes in the Einstein-Maxwell 
theory are characterized by their mass M and charge Q, and exist only for Q < M. 
For Q > M the spacetime describes a naked singularity. One can ask whether it is 
possible to start with a black hole with Q — M, drop in a test particle with charge 
greater than its mass, q > m, and turn the black hole into a naked singularity. 
Wald showed that this cannot happen [2] . In order that a, q > m test particle may 
reach the horizon, it must be thrown in with sufficient kinetic energy that the mass 
of the black hole increases more than its charge. 

In the presence of a positive cosmological constant A, there is a variation of this 
test. Charged black holes in a de Sitter background have three horizons, an inner 
and outer black hole horizon, and a de Sitter horizon. There are thus two extremal 
limits. One corresponds to coinciding inner and outer black hole horizons, and is 
analogous to the extremal Reissner-Nordstrom solution. The other corresponds to 
coinciding outer black hole and de Sitter horizons. In both cases, exceeding these 
extremal cases results in naked singularities. For the first type of extremal limit, 
it was shown in [3] that it is again impossible to destroy the black hole by sending 
in charged test particles. However, we will show in Sec. 2 that for the second 
type of extremal black hole, one can drop in a charged test particle and exceed the 
limiting value. One thus seems to have a potential violation of cosmic censorship 
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when A > 0. 

To analyze this situation further, one needs to go beyond the test particle 
approximation and include the backreaction of the test particle on the geometry. 
Ideally one would hke to have an exact solution describing two black holes coUiding 
in de Sitter space. Remarkably, a solution of exactly this type has recently been 
found [3]. The solution describes an arbitrary number of Q = M black holes with 
A > 0. The solution is dynamical, and the black holes collide in the future. (Alter- 
natively, one can consider the time reverse, which describes white holes splitting 
and expanding.) For small M, the single Q — M black hole is not extremal and 
has three horizons. However, when one increases M, there is an extremal value 
where the outer black hole and de Sitter horizons coincide. One can thus consider 
two black holes each with mass less than this extremal value, but whose sum is 
greater, and let them collide. The purpose of this paper is to examine this process 
and determine whether cosmic censorship is indeed violated. 

In order to settle this question, several preliminary steps are necessary. Some 
of these have interesting applications independent of cosmic censorship. First, one 
needs to investigate the global structure of the multi-black-hole solutions. This 
requires extending beyond the coordinate patch in which the solutions were first 
presented. Here we find a surprise. The horizons across which one must extend 
are in general not smooth, but have only finite differentiability. As a result, the 
extensions are not unique. This has implications for the instability of the inner 
black hole horizon [4, 5], and yields a possible analogy to the lack of smoothness of 
null infinity in asymptotically fiat spacetimes [6]. Physically, for certain horizons, 
the lack of smoothness can be interpreted as the result of electromagnetic and 
gravitational radiation in the spacetime which is not smooth at the horizon. We 
also study null geodesies and determine the global event horizon. We find that 
when the total mass is less than the extremal limit, the black holes coalesce. 

Since A > 0, these solutions are not asymptotically fiat, and the usual def- 
initions of black holes and naked singularities are no longer valid. However, for 
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spacctimcs that asymptotically approach de Sitter, there is an analog of future null 
infinity, and one can define the event horizon to be the boundary of the past of this 
null infinity. Since the singularities of a charged black hole are timelike, they are 
always locally naked. However, if they are enclosed inside an event horizon, there 
is no violation of weak cosmic censorship. We will be interested in whether the 
singularity is visible to all observers who start from a given asymptotic de Sitter 
region in the past. If so, there would be no way for these observers to avoid seeing 
the singularity, and even weak cosmic censorship would be violated. 

Several properties of black holes in de Sitter space have recently been studied 
by Shiromizu et. al. [7]. They discuss black hole collisions, but without the 
benefit of an exact solution. By assuming cosmic censorship, they are lead to the 
conclusion that large black holes will not collide when A > 0. We will see that the 
exact solution behaves quite differently. 

By studying this solution, we find that the question of whether cosmic censor- 
ship is violated by charged black hole collisions is rather subtle. Three different 
issues must be dealt with. First, since the total mass exceeds the extremal limit, 
there is always a naked singularity "on the other side" of de Sitter space, even 
before the black holes collide. Thus, one is essentially starting with singular initial 
conditions and should not be surprised if cosmic censorship is violated. We will 
show that this problem can be alleviated by introducing a charged shell of dust, 
which will remove the unwanted singularity. Second, there is always a Cauchy hori- 
zon, which we will show is generically singular. However, the singularity is rather 
mild, and geodesies can pass through it. Furthermore, there is a large subset of 
the initial data for which the curvature singularity is removed. In these cases, all 
observers cross the Cauchy horizon and see a naked singularity. Cosmic censorship 
is violated. However, the known exact solutions do not describe the most general 
black hole collision with a positive cosmological constant. It is certainly possible 
that in the generic collision the singularity at the Cauchy horizon will be much 
worse, and there will be no way to continue the spacetime. Finally, the naked sin- 
gularity appears only for eternal black holes. We will see that one can form these 
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black holes from regular initial conditions (e.g. using charged dust). But in this 
case, the matter forming the black holes and the matter introduced to remove the 
unwanted singularity collide before a naked singularity is reached. We do not have 
the exact solution beyond this point, and so do not know if naked singularities will 
form in the future. 

Analysis of the global structure of the multi-black-hole solutions is comphcated 
by the fact that, in general, they are neither static, nor spherically symmetric. In 
light of this, we begin in Sec. 2 by discussing the special case of a single black hole. 
This is described by the Reissner-Nordstrom-de Sitter solution. We will show 
that, in the test particle approximation, a Q = M black hole can have its mass 
increased past the extremal limit of coinciding black hole and de Sitter horizons. 
In Sec. 3 we start our discussion of the multi-black-hole solutions by examining 
the basic features of their geometry. By locating the trapped surfaces, we see here 
the first evidence that some of them describe coalescing black holes. In Sec. 4 we 
investigate extensions of the solutions and discuss the lack of smoothness at the 
horizons. Sec. 5 contains a discussion of the event horizon and establishes that 
(when the total mass is less than the extremal limit) its topology indeed changes, 
showing that black holes do combine. In Sec. 6 we consider the introduction of 
charged dust. Finally, in Sec. 7, the possible counterexample to cosmic censorship 
is studied. Sec. 8 contains some concluding remarks. 

2. The Q = M Reissner-Nordstrom-de Sitter Solution 

The Reissner-Nordstrom-de Sitter (RNdS) solution is static and spherically 
symmetric, and its global structure can be analyzed by general methods. This has 
been done by Brill and Hayward [8]. Here we will briefly recount their results for 
the case Q = M. We also discuss briefly the motion of test particles in Q = M 
RNdS backgrounds. We show that one can add q = m test particles to an extremal 
black hole, causing it to exceed the extremal limit. 



5 



2.1. Static Coordinates 



The static form of the Q — M RNdS metric and electromagnetic potential is 

ds^ = -V{R)dT^ + ^ + Fedn\ 
^ ' V{R) 



Wc assume A > 0, and will interpret M as the mass of the black hole. There is 
a curvature singularity at i? = 0. Horizons occur where V{R) = and R > 0. 
The number of horizons depends on the mass and cosmological constant. For 
M^A < jQ, there are three horizons: inner and outer black hole horizons and a de 
Sitter horizon. If the mass of the black hole is increased, the outer black hole and 
de Sitter horizons move closer together. They coincide at the extremal mass given 
by 

ML (2.2) 

For M > Mgxt there is only a single horizon. We will call this the "overmassive" 
case. 

The Penrose diagrams for these three cases are shown in Figs. (1 - 3). Fig. 
(1) shows the case M < M^xt with its three horizons. Notice that are spacelike 
when A > 0. The overmassive case is shown in Fig. (3), and clearly has naked 
singularities. However, these singularities exist for all time, and hence this is not 
a violation of cosmic censorship. The extremal limit. Fig. (2), is rather unusual. 
This solution has singularities, which are visible from $5+, and also a nonsingular 
spacelike surface to their past. (By making an appropriate identification this sur- 
face can even be compact.) However, there are complete timelike geodesies along 
which the singularity is never visible. (These reach the point p in the figure.) Thus, 
one can view the point p as a future timelike infinity and define a horizon to be the 
boundary of its past. This illustrates the difficulty of defining a naked singularity 
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when the cosmological constant is nonzero. As stated earher, we will call a singu- 
larity naked if it is visible to all observers originating from a given asymptotic de 
Sitter region in the past. So the extremal limit is not naked in our sense. 



2.2. Cosmological Coordinates 

The Q = M RNdS solution can also be expressed in isotropic or "cosmological" 
coordinates: 



\jA \ J r U 



where H = ± y^A/S. The transformation from the static coordinates to cosmolog- 
ical coordinates is given by 

1 dh HR^ 

R^Hrr + M, T ^ — log Ht - h(R), ^ = . ^. (2.4) 

H ^ ^ ^' dR {R-M)V{R) ^ ' 

Since (2.3) is the form of the solution that is easily generahzed to several masses, it 
is useful to understand how these cosmological coordinates cover the spacetime. In 
Fig. (1), a single patch of (r, r) coordinates with H < covers the region enclosed 
by the bold lines. Some constant-r surfaces have been drawn in this region. Below 
the dashed line r is negative. The left-hand end of the constant r surfaces is 
r = oo. The right-hand end is r = (but we will soon see that this is an "infinite 
throat," not a regular origin of polar coordinates). Above the dashed line r is 
positive. In this region, the right-hand end of the constant-r surfaces is still r — 0, 
but the chart ends on the left at a finite r, where U vanishes and the geometry is 
singular. Note that this coordinate patch covers a past de Sitter horizon and the 
inner and outer black hole horizons. The future de Sitter horizon corresponds to 
r = oo, T = 0, and the region beyond this horizon, where Q'"'" lies, is not covered 
by these coordinates, but it is covered by (2.3) with H > 0. 

The extremal hmit and overmassive cases. Figs. (2) and (3), are similar. Again 
we have drawn some constant-r surfaces and again r = is shown as a dashed 
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line. The geometry of the spatial slices for r < is regular, but for r > 0, r again 
extends only to a finite maximum value at the singularity. 



2.3. Paths of Test Particles 

Consider a q — m test particle with conserved energy E — m moving in a 
Q = M RNdS spacetime. It was shown in [3] that the radial motion, in static 
coordinates, is given by 



where A is the proper time along the path. There are two possibilities for the 
motion depending on whether the signs in (2.5) are chosen to be the same or 
opposite. The other choices of sign in (2.5) reverse the flow of time along these 
paths. It is simple to check, using (2.4), that choosing both plus signs in (2.5) 
gives a path of constant r in cosmological coordinates. Similarly, we can define a 
new set of cosmological coordinates (r', r') by reversing the sign of the static time 
coordinate T in (2.4). The paths (2.5), with opposite choices of signs are then 
paths of constant r'. These paths of constant r' will be important to us in Sec. 6 
when we discuss the collapse of charged shells. 

Now consider a test particle on a path of constant r, such as the one shown 
in Fig. (1). At early times the test particle is outside both the de Sitter and 
black hole horizons. As time progresses, it enters first the de Sitter and then the 
black hole horizon. Nothing in the above discussion relied on the black hole being 
undermassive. If the black hole is an extremal one, the description is roughly the 
same, with the test particle passing through the degenerate horizon. It follows 
immediately that there is no barrier to such a test particle entering the horizon. 
One can exceed the extremal mass limit by dropping in test particles. Furthermore, 
since the charges on the black hole and test particle provide a repulsive force, we 
expect that this phenomenon will be generic for black holes and test particles with 
charge less than their mass. 





(2.5) 
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3. General Properties of the Multi-Black-Hole Solutions 



3.1. Geometry of spatial surfaces 

The solutions of reference [3] depend on a number of parameters that corre- 
spond to several different masses at arbitrary positions (but not arbitrary veloc- 
ities). The metric and gauge field for mass parameters Mi and positions fi are 
given by 

ds^^-^ + U'^df ■ df, Ar^\. 

M- (3-1) 



I 



Unless otherwise stated, we will assume H < 0. It is only for this case that the 
spacetime describes black hole collisions. In terms of H, the extremal limit (2.2) is 

Mext = (3.2) 

The surfaces of constant r are spacelike everywhere. Near each fi the geom- 
etry resembles the infinite throat familiar from the asymptotically flat extremal 
Reissner-Nordstrom solution. This can be seen by expressing the spatial metric in 
spherical coordinates centered at r = r^. Near the origin of these coordinates, this 
metric becomes 

Mf 
r2 



df ^ + Mfdn^, (3.3) 



which is the metric for a cylinder of infinite spatial extent having cross sectional 
area AnMf. 

The curvature of (3.1) can be singular at zeros of the metric function U. This 
can be seen from the square of the Maxwell field strength, 

f2 = F^.F^'^ = (3.4) 

— * — * 

where V denotes the usual gradient with respect to f. li U — and VC/ does not 
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vanish like f/^ or faster, then diverges and the curvature is singular. The metric 
is regular everywhere else, but it is incomplete as discussed in Sec. 4. 

Assembling these elements we get a picture of how the spatial geometry de- 
velops in cosmological time. We will start with one black hole first. Consider the 
spatial surfaces for a single black hole with mass M < Mg^t, which are sketched on 
Fig. (1). The metric function is U = Hr + For r < 0, f/ is positive everywhere 
and the spatial surfaces are nonsingular. They are asymptotically flat and have the 
cylindrical form of an infinite throat near the origin. For r — 0, the spatial metric 
is regular and has the cylindrical form everywhere. As r is increased slightly, a 
singularity appears near r — oo. As r increases further, the singularity cuts off 
more and more of the cylinder. In Fig. (1) this is shown by the way the singularity 
intersects the spatial surfaces. 

The generalization to more than one mass is then straightforward. For r < 
the spatial surfaces are nonsingular and asymptotically flat at large radius. Near 
each "point" rj, the spatial metric has the form (3.3) of a throat. These surfaces 
are depicted in Fig. (4a). For r — 0, the surface is nonsingular, but spatial infinity 
is now also asymptotically cylindrical. This surface is depicted in Fig. (4b). As 
r is increased slightly above zero, a singularity moves in from spatial infinity as 
shown in Fig. (4c). As r continues to increase, the singularity splits and eventually 
surrounds each of the throats individually. The spatial surface is then composed 
of a number of isolated throats, as depicted in Fig. (4d). This description is 
independent of the size of the masses, and hence applies to both the overmassive 
and undermassive cases. 

3.2. Trapped Surfaces 

The causal structure of the multi-black-hole solutions is rather involved. To 
gain some understanding, we begin by considering trapped surfaces. The expan- 
sions Bout {din) oi outgoing (ingoing) null rays normal to a 2-surface E in a spatial 
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hypersurface are given by 



9 = Dy - Kiju'n^ + K, (3.5) 

where hij, Kij are the metric and extrinsic curvature on the spatial hypersurface, 
Di is the covariant derivative compatible with hij, K — hijK'^^ , and is the 
outward (inward) directed normal vector to E. A surface E is called outer trapped 

if 6out < 0, and inner trapped if > 0. Surfaces for which Oout (Oin) vanish are 
called outer (inner) apparent horizons. For the surfaces of constant r in the metric 
(3.1), the extrinsic curvature is simply given by Kij = Hhij. 

In the case of a single black hole of mass M, one can calculate the expansions 
exactly for spheres centered on the origin. The result depends only on the quantity 

R = Hrr, (3.6) 

which is simply related to the static R coordinate by R — R — M (2.4). The 
expansions are 

It will be useful to define 

a = Vl + ^MH (3.8) 



and 



dout vanishes at 



and 



/3 = Vl - 4:MH. (3.9) 



Rbh^-^{l + 2MH-a) (3.10) 



RdS^-^{l + '2MH + a), (3.11) 
which correspond to the black hole and de Sitter horizons. The ingoing expansion 
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6in vanishes at 

Rin = ^{l-2MH-P), (3.12) 

which corresponds to the inner horizon. The quantities a and (3 will play an 
important role in what follows. 

In Fig. (5) we have drawn the coordinate patch covered by (r, r). The co- 
ordinate patch is divided into four regions, labeled I - IV, by various horizons. 
These horizons coincide with the boundaries of the regions of trapped surfaces, 

Rds, Rim given in (3.10) - (3.12). Spheres in region I are outer trapped, and 
the boundary between regions I and II corresponds to Rds- Spheres in region III 
are also outer trapped and the boundary between regions II and III is given by 

Spheres in region IV are both inner and outer trapped. 

We now consider the two-black-hole solution where each mass is less than the 
extremal limit Mi < M^xt- The main difference from the one-black-hole case, or 
the multiple black holes with zero cosmological constant [9] , is that these solutions 
are dynamical. The apparent horizon will evolve in time. Of course, with two 
black holes an apparent horizon will no longer be precisely spherical, but it will be 
approximately spherical in limiting cases. Hence we consider 

eout = 0, (3.13) 

evaluated on appropriately centered spheres. (Trapped surfaces have been found 
for initial data describing two uncharged black holes with A > in [10].) 

Spheres of sufficiently small radius, centered at fj, are always trapped. At early 
time, for r <^ 0, there is a solution for 6out = at |r — fi\ = Ri,j^[Mi]/ Ht, since 
the equation is identical to that of RNdS with mass Mj. (Here i?ft/j[Mj] is defined 
as in (3.10), with M replaced by Mj.) So, around each throat there is a region of 
outer trapped surfaces, surrounded by an external region. Actually, around each 
mass there is also the second solution to 9out = at |r — fi\ — Rds[Mi\/ Ht. These 
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two solutions correspond to the black hole horizon and the past de Sitter horizon 
in the RNdS spacetime. 

Spheres of constant r (where r is much larger than the coordinate distance 
between the two masses) resemble spheres in the RNdS solution with mass equal to 
the total mass M — Y^Mi. If M < M^xt, then at late times, as r approaches zero, 
one again has two (outer) apparent horizons at r = Rbh/Hr and r = Rds/Hr. This 
suggests that the black holes coalesce. If M > Mg^t: then there are no apparent 
horizons at late time for large r. We will have more to say about this case in Sec. 
7. 

One can understand the behavior of the apparent horizon as follows. For the 
RNdS solution, the horizon is at fixed R, which corresponds to r oc 1/Ht. So when 
T is large and negative, the horizon is at small r, and when r ~ 0, it is at large 
r. In both of these limits, the two black hole solution resembles the one black hole 
solution (although with different masses). 



4. Extensions 

4.1. Locating the Horizons 

As one might expect from the RNdS solution, the region of spacetime described 
by the metric (3.1) (with r, f taking all real values) is incomplete, even away from 
the curvature singularity. It is bounded by the analog of the de Sitter horizon at 
large r and the inner black hole horizon and past white hole horizon at small \f—fi \ . 
In both these regions the metric becomes approximately spherically symmetric. 
We will establish the incompleteness by considering radial null geodesies in these 
asymptotic regions. 

We first derive an equation for the affine parameter s along a radial null geodesic 
from the variational principle 
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Vary r and use the outgoing null condition, 



dr/dr = U^, (4.2) 

to find 

2HU(^X = 0. (4.3) 



If a null geodesic is known as r = T(r), the affine parameter is given by quadratures, 

j e^HiU{r(r),r)dr^^ (44) 

In the limit of large r, the function U takes the simple form 

M 

U^Hr + —. (4.5) 
r 

where M — '^ Mj. In this case, the null condition (4.2) is integrable when rewritten 
in terms of the variable R of (3.6): 

^ =R + H(R + Mf. (4.6) 
dr 

We need only the asymptotic form of the solution. If M < Mgxt, then the right 
hand side has two roots at Ri^h and Rds- Starting with R > Ri,i^ one finds that the 
solution always approaches Rds as r — >• cxd (see Fig. 6). Thus 

Rds + M 

U {T{r),r) —>■ . (4.7) 

r 

Using this in (4.4), and noting from (3.11) that R^s + M = —(1 + a)/2H we find 
that 

sr^ SH + cr'"" (4.8) 

where sh is the horizon value of s, and c is a constant of integration. Thus s 
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remains finite as r — oo. Asymptotically for large r, it follows that 

r- (s-s//)-^/", r-(s-Siy)V". (4.9) 

This not only shows that these quantities reach their horizon values at a finite s, 
but also that 1/r and r are not smooth functions of s (unless 1/a happens to be 
integral) . 

For each mass Mi, we can find similar results for ingoing null geodesies near the 
inner horizon, in coordinates centered about Mj. In the limit of small r, ingoing 
null geodesies satisfy 

The last term on the right is a constant which can be removed by shifting the 
origin of r. We then obtain an equation like (4.6) whose solution has the hmiting 
form 

R^Rinm, £/(r(r),r)^ ^-^^;^ + ^ . (4.11) 
Eq. (4.4) now implies that as r goes to zero, 

sr^si + Cr^' (4.12) 

where sj denotes the value of the affine parameter at the inner horizon, C is a 
constant of integration and Pi = y^l + 4Mj|iJ|. So the affine parameter is again 
finite at the inner horizon. It follows that 

r-(s-S7)^/^% {s-si)-^/^\ (4.13) 
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4.2. Extending Beyond the Horizons 



We now consider extensions across the horizons. We begin by introducing 
new variables. For the single black hole, the static R coordinate is good on the 
horizons. The closest analog for the solutions (3.1) is R = Hrr introduced above, 
where r = is chosen to correspond to the location of one of the masses. We also 
set 

y = lnr, W = rU. (4-14) 
In {R,y) coordinates, the metric (3.1) takes the form 



- - ^'^^^^f + WW + (4.15) 



We have seen that the metric approaches the solution for a single black hole 
both in the limit of large and small r. However, in determining the behavior of the 
geometry across the horizon, the rate at which the metric approaches the single- 
black-hole solution is crucial. We first show that all curvature scalars remain finite 
as one approaches these horizons. The (i?, y) part of the metric (4.15) has constant 
determinant. Thus the inverse metric has a similar form, with only W'^ appearing 
in the denominator. A general curvature scalar will involve terms consisting of 
derivatives of the metric and its inverse, multiplied by powers of the metric and its 
inverse. All of these terms reduce to derivatives of W and R divided by powers of 
W. But all derivatives of W remain bounded as y ^ ±oo, and since W is finite 
on the horizon, these terms cannot blow up. (By simply shifting the origin of r, 
one can apply this argument to the horizons near each of the masses.) 

It is tempting to conclude from this that the metric is smooth across the hori- 
zon and can be analytically continued as in the single-black-hole case. However 
this is incorrect. We will see that, in general, the horizon has only finite differen- 
tiability. The curvature can even be singular at the inner black hole horizon, but 
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the singularity is "null," and so all curvature scalars remain finite . To establish 
this result, we will introduce coordinates that are good in a neighborhood of the 
horizon. It will turn out that i? is a good coordinate but that r = e^^ is not. (We 
have already seen in (4.9) and (4.13) that r is not a smooth function of the affine 
parameter along null geodesies.) For the single-mass solution, W — R + M, so the 
metric (4.15) depends only on R, and therefore is smooth. However, the effect of 
the other masses is to modify the single- black- hole metric by a power series in r, 
so the exact metric is no longer smooth across the horizon. 

The metric (4.15) can be rewritten in the form 

ds'^ = + ^^^^ ~ ^^^^] + ^^^^ + + W^dn'^. (4.16) 

Our procedure for studying the extensions of the spacetime is to introduce new 
coordinates {u, v) in the neighborhood of each horizon, such that the first quantity 
in brackets is proportional to du and the second is proportional to dv at the horizon. 

We first consider the de Sitter horizon {y — > oo). If we choose the origin of 
coordinates to be the center of mass, the metric takes the form (3.1) with 



U^Hr + ^ + l^ + 0{r-') (4.17) 



where M is the sum of the individual masses. In terms of {R, y) coordinates the 
metric takes the form (4.16) with 

W ^ R + M + f{e,(t))e-'^y + (4.18) 

We saw in Sec. 4.1 that R approaches the value R^is at the horizon. Expanding 



★ The above argument shows that aU components of the Riemann tensor in (i?, y) coordinates 
stay bounded, but this coordinate basis is not well behaved at the horizon. 
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HW"^ about Rds yields 



HW^ ^ -Rds -il + a)[R-Rds + /e^^y]. 



(4.19) 



Define new coordinates 



u= 1 + 



{R - Rds) + 



a{a + l)f 



-2y ^-ay 



(4.20) 



'2Rds 



2Rds{a + 2) 



v^- R-Rds + 



{a + l)f 



a-2 



The metric near the horizon becomes 



2Rds 



dudv + W^dn^. 



(4.21) 



So u, V are good null coordinates near the horizon, which now corresponds to 
u = 0. The curves v = constant cross the horizon. The metric (4.21) depends on 
W'^ which involves factors of (R — Rds) and e~^^. In addition, there are corrections 
to the leading order behavior (4.21) that involve these same factors. To express 
these terms as functions of u and v we need to invert (4.20). Near the horizon 
we have u = e""?^, and R - Rds ^ -[ve-^f + {a + l)fe-'^v/{a - 2)]. But since 
a — sjl — AM\H\ < 1, as long as 7^ the second term is neghgible compared 
to the first as y — > 00, and we have R — Rds — uv. Thus corrections involving 
only powers oi R — Rds will be smooth at the horizon. However, e" -y = u^/^ is 
not smooth in general. Since the corrections start at order e~^^ = u^/** and a <1, 
we see that the metric is always at least in these coordinates. Since the metric 
is C^, there is no curvature singularity. But one can show that in general, certain 
components of derivatives of the curvature diverge at the horizon, so there are no 
coordinates for which the metric is analytic. As a result, the extension across the 
horizon is not unique. One can match onto essentially any solution of the form 
(3.1) with the same total mass. (Similar behavior of finite differentiability across 
a horizon in an exact solution was found in [11].) 
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The differentiability at the horizon can be increased in two ways. If the total 
mass M is chosen so that 1/a is an integer n, then e~y is smooth and the metric 
is C°° . This occurs when 



For these values, the smooth continuation consists of matching the spacetime onto 
one with the same positions and magnitudes of all the masses so that all multipole 
moments agree, but with the opposite sign of H. We do not understand the physical 
significance of these special masses. 

The second way to increase differentiability is to arrange the masses so that 
their first n multipole moments vanish. Then the perturbation will begin at e^""^. 
These solutions may provide a simple model of smoothness of null infinity in asymp- 
totically flat spacetimes. It has been suggested [12] that the behavior of fields at 
null infinity may depend on their fall-off near spatial infinity. If the unconstrained 
part of the initial data falls off more quickly in spacelike directions, then perhaps 
the evolved fields will be more differentiable at null infinity. This is very similar to 
the behavior we find in the multi-black-hole solutions. 

We now consider the inner horizon associated with one mass Mj which we 
choose to be located at the origin. Near r = we have 



where c is a constant that can be removed by changing the origin of r. Thus 
W = R + Mi + gr'^ + ■ ■ ■. We showed in Sec. 4.1 that R approaches the constant 
Rin[Mi\ at the inner horizon. Expanding HW^ near R = Rin[Mi\ yields: 



4.M\H\ = 1 - — . 



n- 



(4.22) 



U = HT+^ + c + g{e,(P)r + ■■■ 



{A.2'i) 



HW^ ^ Rin[Mi] + (1 - A) {R - Rin[Mi] + ge^y) 



(4.24) 



where, as before, Pi = ^1 + AMi\H\. Define new coordinates by 
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Then the leading order behavior of the metric near the horizon takes the simple 
form 

(The leading factor in parentheses is just a constant.) So u,v are again good 
null coordinates near the horizon, which now corresponds to f = 0. Near the 
horizon we have v = e^'^, and R — Rin = ue^'^ + {Pi — l)ge^y / (2 — j3i). But since 
every black hole with Mi < Mext has (3i < 2. Thus as long 
as 7^ the second term is negligible compared to the first near the horizon, and 
we have it! — Rin[Mi\ — uv. Corrections involving only powers oi R — Rin[Mi\ will 
again be smooth at the horizon. However, r — — and 1 < A < 2. Thus 

r is not smooth at the inner horizon. Since the corrections start at order r^, the 
metric is but not in these coordinates. 

To see that there are no better coordinates for which the metric is smooth, 
we compute a component of the curvature. Let I — d/dv and 77 = d/dO. Then 
RyOve — Riivpol^'H^l^'n'^ contains several terms that are finite at the horizon. But 
it contains one term which is infinite there. This is gee,vvi which involves two 
derivatives of v^/^\ This divergence is null since it only occurs at v = 0. We will 
return to the physical interpretation of this singularity shortly. 

Unlike the case of the de Sitter horizon, there are no special values of the mass 
for which the inner horizon becomes smooth. However one can still increase the 
differentiability by carefully arranging the other masses so that the first n powers 
of r cancel in the expansion of U in (4.23). In particular, one could remove the 
curvature singularity this way. However, the inner horizon associated with the 
other black holes will still be only C^. 

It is interesting to compare this situation with the case of zero cosmological 
constant. There, it was shown [13] that the spacetime describing several Q — 
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M black holes has an analytic extension across each inner horizon, which simply 
corresponds to letting r become negative. In other words, when A = 0, r is smooth 
at the horizon and there is no analog of the singularity we find above. This is 
consistent with our results since A = implies f3i — 1 for all i, so r = at the 
inner horizon. 

The remaining horizon is the past white hole horizon r = 0, t — — oo. One 
can show that the situation here is similar to the de Sitter horizon. The metric is 

always and can be made C°° if the individual mass Mj takes one of the discrete 
values (4.22) . Notice that if more than one mass has one of the discrete values, 
the total mass exceeds the extremal limit (3.2). 

4.3. Physical Interpretation of the Lack of Smoothness 

We now consider the physical interpretation of the lack of smoothness we have 
found at the horizons. Consider the solution with two black holes and take the Hmit 
where one mass becomes much less than the other. In this limit, we can think of the 
small mass as a test particle moving in the background of the large black hole. The 
single black hole solution is, of course, smooth everywhere. But the perturbation 
in the metric and Maxwell field obtained in this limit involves a function of r 
(or 1/r), which is not smooth at the horizon. This seems rather unphysical. We 
would expect a q — m test particle to radiate as it falls into the black hole, and 
the field it produces should remain smooth at the de Sitter horizon. This can be 
made more precise in terms of initial data. To avoid singularities associated with 
point particles, we can model the test particle by a small ball of g = m dust. Let 
us take our initial surface to be the T = surface in the RNdS metric. One can 
certainly find smooth initial data for the linearized Einstein- Maxwell field equations 
that satisfy the constraints with the ball of dust as a source. The evolution of 
these initial data must be smooth everywhere in the domain of dependence, which 
includes both the de Sitter horizon and past white hole horizon. The fact that the 
multi-black-hole solutions are not smooth at these horizons can be interpreted as 
saying that they describe more than colliding black holes. In addition, they contain 
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a distribution of electromagnetic and gravitational radiation which is not smooth 
everywhere in the spacetime. We expect that there are other (undoubtedly more 
complicated) solutions describing colliding black holes that do not suffer from this 
lack of smoothness. 

The singularity at the inner horizon is qualitatively different. This is because it 
is on the boundary of the domain of dependence of the initial data surface described 
above. The behavior we find is reminiscent of the instability of the inner black hole 
horizon which has been extensively studied for A = [4]. (When A > it has been 
argued that the inner horizon might be stable [5].) A key difference, however, 
is that the previous analyses were based on a perturbation expansion, while the 
solution (3.1) treats the effect of the other masses exactly. One unusual feature 
of these solutions is that only half of the Cauchy horizon becomes singular. The 
other half lies inside the (r, r) coordinates and is C°°. Perhaps this is related to 
the extra radiation in the spacetime that is responsible for the lack of smoothness 
at the de Sitter horizon. 

The idea that the finite differentiability at the de Sitter and past white hole 
horizons is a result of additional radiation is supported by the fact that one can 
construct exact initial data for multiple charged black holes with Q < M. These 
initial data are smooth at the white hole horizon. This construction relies on 
the A = initial value solutions of Brill and Lindquist [14] for time-symmetric, 
arbitrarily placed wormholes of general mass and charge, and on the method of 
Nakao et al [10] to turn such a solution into one with a cosmological constant. The 
3-metric and electromagnetic potential of [14] are 



where u is a, unit normal form to the initial surface, and a, h and c determine the 
wormholes' masses and charges (approximately, Mj ai/c + bi, Qi ~ —ai/c + bi). 



ds"^ = {x^pfidx^ + dy'^ + dz^) A = u ln(x/^) 



(4.27) 
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Eq (4.27) satisfies the A = Einstein constraints if the extrinsic curvature Kij 
vanishes; therefore it satisfies the A 7^ constraints if we put [10] 



Kij = ±Hgij, H = 7X73. (4.28) 



Special cases of (4.27) are initial values of the known dynamical solutions, as 
follows: The general Q ^ M RNdS geometry in cosmological coordinates [8], 



t/2 

(1 - ^) (dR - Rdyf 



+ W'^{dy'^ + d^l^) (4.29) 



U -I ^ -Q^ W- R M 

has the form (4.27) with ai = (M - Q)/2, 61 = (M + Q)/2, = = 6j for i ^ 1, 
c = 1, on the surface t — Q. The solutions (3.1) have the form (4.27) with = 0, 
hi — mi, and c = Ht. 

We also note that the geometry of the surface r = of the solution (3.1), as 
well as that of the surface R = M ± a/ — Q2 ^j^g general RNdS geometry 
(in static coordinates) can be generalized by initial values of the type (4.27) with 
c = 0. 
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5. Event Horizons 



All of the extensions considered above have a S^"*" beyond the "de Sitter hori- 
zon," which corresponds to the limit (r, r) (oo,0), but R — > R^s finite, in our 
original coordinates (3.1). In these original coordinates we can therefore identify 
null curves that can go to as those that reach this horizon. By contrast, the 
black hole horizon as defined in the Introduction (see also Ref. 7) is contained 
entirely within the original chart: it is the boundary of the region of events that 
can be causally connected to the de Sitter horizon. Thus we identify events from 
which outgoing null geodesies must reach the U = singularity as lying inside the 
black hole horizon. 

In this section we establish several properties of null geodesies that have im- 
plications for the black hole horizon. For simplicity we confine attention to the 
solution (3.1) for two centers with identical mass parameter M/2 with M < M^xt, 
separated in the 3D Euclidean base space by distance 2a. The origin of the Eu- 
clidean coordinate system is at the midpoint between the two centers. The Eu- 
clidean line connecting the two will be called the axis, and the perpendicular plane 
at the origin is the midplane. We show that for each sufficiently late (but negative) 
time r there is a sphere of radius r(T) such that all outgoing null geodesies from 
the sphere will reach the U = singularity — hence no causal curve from inside 
can reach r = oo; but that there arc points outside the sphere from which causal 
curves reach r = cxd in finite affinc parameter. That is, the sphere lies within the 
event horizon. We also show that at sufficiently early times, all points in the mid- 
plane can be connected to r = oo by causal curves; in other words, at early times 
the horizon does not meet the midplane. Finally, we will show that at sufficiently 
early times the event horizon consists of a sphere centered on each mass. Taken 
together these results prove that the black holes coalesce. The event horizon at 
early times has two disconnected pieces, but at late times it has only one. 

We first consider the region inside the horizon. Every outgoing null geodesic 
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must satisfy 



UV < ^. (5.1) 



A lower limit on the "potential" U is obtained by pretending that the total mass 
M is concentrated at the greatest possible distance, r + a, so that 

^>C2>fHr + ^y. (5.2, 

In terms of new variables Rk — Hrir + a), = ln(r + a) the inequality simplifies, 

JIT3 

-^<R^ + H{IU + Mf = H{R^ - Rds)(R* - Rbh)- (5.3) 

The slope of the function on the RHS is positive at its smaller zero, R^h (Fig. 6). 
Suppose that initially il^ < R^h, i-e., a < r < (Rhh/Hr) — a. This can be fulfilled 
if T > —R}}ii/2a\H\. Then R^ decreases with r and becomes negative at a finite r. 
But then r is positive, and as it increases further with r, the singularity at C/ = is 
reached at a finite r. Thus, as r approaches zero (from below), there exist spheres 
of constant r enclosing both masses such that all outgoing null geodesies hit the 
singularity. 

We now show that at early times, a similar result holds for small spheres 
centered on each mass. For simplicity, let us shift the origin of spherical coordinates 
to be at one of the masses, and assume r < a. Then the "potential" from the other 
mass is at least M/6a, so that instead of Eq (5.2) we have 



-r>U^> 
dr 



^'^+6^) + 27 



(5.4) 



Introducing the new variable R+ — Hr{r + g^) leads to an equation analogous 
to (5.3), 

dR+/dy <R+ + H{R+ + \Mf. (5.5) 

Analyzing (5.5) in the same way as (5.3) we then find that events satisfying Rj^ < 
Rl,h[M/2\ are inside the horizon. This inequality is satisfied at any early (negative) 
T by all sufficiently small r < Rbh[M/2]/{HT + M/6a). 
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To explore the outside of the black hole horizon, we ask under which conditions 
some causal curves from a given event can reach r = oo. We first confine attention 
to the midplane, and obtain an upper limit on the inverse speed of purely radial 
null curves: 

= U^^ lHr + -^=== <\Ht + . 5.6 

dr V Vr^ + aV V + 

In terms of the variable i?^, introduced in Eq (5.3) this becomes 

dR^/dy^> R^ + H{R^ + V2Mf. (5.7) 

Arguing as before (see Fig. 6) we find that if i?^ > Ri,h[V^M] is satisfied initially 
for some r, then will stay positive for all larger r. But this means that r will 
remain negative as r increases, avoiding the singularity. 

For example, if the null curve starts at any r, but with r < Rbh[V2M]/Ha, 
it satisfies the initial inequality and hence will be able to escape to r = cxo. Thus 
at early times no point in the midplane is within the horizon (provided of course 
that a real Rf,h[V2M] exists; if M is close to the extremal hmit, a more accurate 
investigation of the null geodesies is needed) . 

Similarly, for any r < we can find a sufficiently large 

r> {Ri,h[V2M]/HT)-a (5.8) 

so that Ri, > Ri,fi[V2M] and null geodesies can escape to infinity. In fact, for 
sufficiently large r the argument leading to the increase of i?* is valid also for points 
not on the midplane. Thus outside the sphere of events that must causally lead 
to the singularity there are events from which causal escape to infinity is possible. 
At early times, therefore, the event horizon surrounds each mass separately, and 
expands approximately according to r = Ri,h[M/2]/HT. At r ^ Rijh/2aH the 
event horizon enters the midplane, and thereby changes its topology to a surface 
surrounding both centers. 
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6. Charged Shells and Dust 



As stated in the Introduction, in our test of cosmic censorship, we will need to 
introduce a shell of charged dust to remove an unwanted singularity. The dynamics 
of these shells is also of interest in its own right. For example, one can ask whether 
it is even possible to form Q = M RNdS black holes from collapsing matter. For 
A = 0, Boulware [15] studied the dynamics of charged shells. He found that a 
shell having charge density equal to its rest mass density and also equal to its total 
mass density does not collapse. Rather, as one might expect, it stays at constant 
area. For A > 0, on the other hand, we find that a simple extension of Boulware's 
calculation does give collapse to form Q — M RNdS black holes. 

6.1. Shells 

Consider a spherical charged shell in an otherwise empty spacetime with A > 0. 
Inside the shell Birkhoff 's theorem guarantees that the metric is de Sitter. Outside 
the shell, we match the metric to RNdS. In order to determine the motion of the 
shell, one integrates Einstein's equation across the shell to obtain jump conditions 
on the curvature. Working in static coordinates, one finds that a shell having 
charge Q equal to its rest mass M also equal to its total mass {i.e. no kinetic 
energy) follows the same path (2.5) as the radially moving q = m test particle 
discussed in Sec. 2.3. In cosmological coordinates then, the shells stay at constant 
comoving radius. 

The physical picture of shell collapse is then quite simple. The metric has 
the form (3.1), with U = Ht ior r < Vg and U = Ht + M/r for r > r^. The 
comoving radius of the horizon however, changes with time. For M < Mg^t, there 
is a horizon at i? = Htjjt = Ri,j^. For large negative r, rjj is inside the shell and 
hence does not correspond to a horizon. As r increases towards zero, ru sweeps 
past the shell, which is now contained within the black hole. The same process can 
be used for a number olQ — M shells to form the multi-black-hole solutions, if the 
shells are located at an equipotential of U and have a mass and charge distribution 
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appropriate to a constant U in the interior. Then the metric is de Sitter inside 
each sheU, and the muhi-black-hole metric outside. AU the sheUs stay at constant 
comoving radius, and eventuaUy become black holes which later merge. Spacetime 
diagrams for a single collapsing shell are given in Figs. (7a,b) for undermassive 
and overmassive shells. The shells take a finite cosmological time to reach the 
singularities, but an infinite proper time. 

The shells in Fig. (7) all start on the right side of the Penrose diagram and 
collapse to the singularity on the left. We can also introduce shells on the left which 
collapse to the singularity on the right. This follows from the T — > — T symmetry 
of the RNdS solution in static coordinates. The motion of these shells is simple in 
the (r', r') coordinates defined in section 2.3. The shells stay at constant r'. For 
r' < r'g the metric will be de Sitter, and for r' > r'g the metric will be RNdS. These 
shells will be useful for us in Sec. 7 in constructing our potential counterexample 
to cosmic censorship. 

6.2. DUST 

It turns out to be easy to handle arbitrary configurations of g = m dust as 
well. Consider a general metric and gauge field of the form 

where dU/dr = H. For matter consider charged, comoving dust. The matter 
stress-energy and current density are given by 

J'dust _ p^^g^UaUf), Ja — Pcharge'^ai (6-2) 

where u"" — jj{-^Y is the four velocity of the dust. The Hamiltonian constraint is 
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given by 

= -lQ'K{pdust + P Maxwell + Pcosmo)- 

The third term on the second hne is just the cosmological term. The second is 
equal to the energy density of the Maxwell field. Thus, the first term must be 
compensated by the energy density of the dust. For a solution then, one requires 



V^C/ = -ATTpdustU^- (6.4) 

The Maxwell constraint gives the same equation (6.4), if Pcharge — Pdust- The 
evolution equations require that T^-^^^ — and that dr{U^Pdust) — 0. The first 
requirement is satisfied for comoving dust. The second condition requires that the 
function p{x) = U^Pdust be independent of r. The constraint equation (6.4) is then 

V^t/ = -Anp. (6.5) 

We see that we can freely specify the function p. This, together with dU /dr — 
H, and the boundary condition U Hr as r ^ oo determine a unique U . The 

physical dust density Pdust is then derived from p and U . Note that the volume 
element \/h = for the spatial metric, so that the total mass of a dust cloud is 
constant and given by 

M = y d^xVhpdust = J d^xp. (6.6) 

Other measures of the dust cloud, e.g. its volume, are of course time dependent 
and consistent with the picture of a collapsing cloud. 
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7. Testing Cosmic Censorship 



Having developed the necessary properties of the multi-black-hole solution we 
now turn to our test of cosmic censorship. Consider two black holes, each with 
mass less than the extremal limit, but whose sum is greater. Since the multi-black- 
holc solutions are not spherically symmetric, it is difficult to describe their global 
structure completely. It is therefore convenient to focus on a two dimensional slice 
of the spacetime. Consider first a curve which comes up one throat and down the 
second on a constant r surface. Evolving in time, we obtain a two dimensional 
subset of the spacetime which is nonsingular for r < 0, but develops a timelike 
singularity between the two throats at a certain time tq > (when U — 0). This 
appears to be a violation of cosmic censorship in which a naked singularity evolves 
from regular initial conditions. However, since this subset of the spacetime does 
not include the asymptotic region, one does not know if the singularity is hidden 
behind an event horizon. 

A better choice is to consider the curve which starts at infinity and goes down 
one of the throats on a constant r surface. Evolving in time, we obtain a two 
dimensional slice of the spacetime whose Penrose diagram is shown in Fig. (8a). 
For small r, the solution looks like a single subcritical black hole, except that 
has been pushed ofi^ to the future of the region shown. For large r, it resembles 
the overmassive case with its singularity. This singularity is naked and exists for 
all time. It is independent of whether the two black holes have collided or not, 
and is just a refiection of the fact that the total mass always was greater than the 
extremal limit. Fortunately, this singularity can be removed by adding a shell of 
charged matter. This follows from the fact that the solution, at large r, reduces to 
the overmassive RNdS solution, and we saw in Sec. 6.1 that the singularity in this 
solution could be removed by a shell. The Penrose diagram for the solution with 
the shell added in shown in Fig. (8b). 

In the absence of a cosmological constant, one usually requires that a coun- 
terexample to cosmic censorship have nonsingular data on a surface that is asymp- 
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totically flat outside of a compact set. In the presence of a cosmological constant, 
the analog would be nonsingular data on a compact manifold. Our example does 
not have a compact surface, but it has what might be considered "the next best 
thing." Consider the surface S shown in Fig. (8b). This surface is defined by r 
equal to a negative constant outside the shell, and any spacelike surface inside the 
shell that continuously joins to it. The initial data on this surface is nonsingular 
everywhere. The surface is not compact, but has two infinite throats (only one of 
which is shown on the figure). However, each of these throats is surrounded by a 
trapped surface. So one would not expect that the asymptotic regions down the 
throats could infiuence the solution in the interior. 

The initial data on S uniquely determines the solution up to the Cauchy hori- 
zon. As we have discussed in Sec. 4, the solution past the horizon is not unique 
but all extensions have a curvature singularity shown at the right in Fig. (8b). It 
is clear from the diagram that all observers originating from 9~ reach the Cauchy 
horizon, and if they extend beyond, they will see the singularity. We have also seen 
in Sec. 4 that the horizon is not smooth. Generically, for the class of solutions 
(3.1), the horizon is but not C^. If we end the spacetime at the horizon then, 
of course, there is no violation of cosmic censorship. But one can find choices of 
parameters such that the Cauchy horizon associated with one of the masses is at 
least and there is no curvature singularity. Furthermore, even in the general 
solution, geodesies can be extended beyond the horizon. So it seems reasonable to 
conclude that cosmic censorship is violated in these examples. 

It is clear that this violation is associated with the infinite throats in the initial 
data. One way to see this is that in cases when the horizon is or smoother, it is 
homogeneous. There is no point where the curvature is becoming large, associated 
with the beginning of a naked singularity. The singularity seems to "come in from 
infinity." The fact that the asymptotic regions are hidden behind trapped surfaces 
in the initial data, does not seem to be sufficient to prevent the violation of cosmic 
censorship. In a sense, all of space collapses down the throats carrying all observers 
with it. 
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In light of this, it is natural to ask whether cosmic censorship would be violated 
if one first formed the black holes from regular initial conditions. We have shown 
that one can, in fact, form the black holes using shells of charged dust. If one also 
removes the singularity at infinity with another shell as we have discussed above, 
one has compact initial data. The problem now is that the two types of shells 
collide in the spacetime before any singularities have formed. We do not know 
the solution explicitly after this occurs. However, even if naked singularities were 
found later in the evolution, one would not know whether they were fundamental, 
or artifacts of the dust approximation. It is well known that naked singularities can 
form in spherically symmetric dust collapse. These shell crossing or shell focusing 
singularities can also occur in the absence of gravity and hence have nothing to do 
with cosmic censorship. 

8. Conclusions 

Motivated by a new test of cosmic censorship, we have studied the global 
structure of the multi-black-hole solutions (3.1). We have found that, if the total 
mass is less than the extremal hmit (3.2), then they describe black holes which 
coalesce. It is remarkable that an analytic solution describing coalescing black 
holes can be expressed in such a simple form. Somewhat surprisingly, we have also 
found that these solutions contain radiation which is not smooth at the de Sitter 
(and past white hole) horizon. Perhaps the presence of this radiation is related to 
the simplicity of the solution. 

The test of cosmic censorship was based on the fact that there is an upper limit 
to the mass of a black hole when the cosmological constant is positive. We have 
seen that colliding two black holes which are each less than the extremal mass, but 
whose sum is greater, does produce naked singularities. 

However, we cannot yet claim that this is a serious violation of cosmic censor- 
ship for two reasons. The first concerns how generic the violation is. As we have 
seen, the most general of the exact solutions has the naked singularity protected 



32 



by a Cauchy horizon with a weak singularity. However, the exact sohitions only 
describe a subset of black hole collisions with a positive cosmological constant. The 
initial position and masses can be specified arbitrarily but not their initial veloci- 
ties. And, of course, one cannot specify arbitrary gravitational and electromagnetic 
radiation. It is not clear whether the most general solution has a Cauchy horizon 
with a stronger singularity. If so, then cosmic censorship would be preserved. 

The other reason concerns the fact that the example involves eternal black holes 
and not ones which formed from compact initial conditions. It is not yet clear how 
physically reasonable collapse would affect the formation of naked singularities. 
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Figure Captions 

Fig. 1: Penrose conformal diagram for the Q — M RNdS geometry with M < 
Mext- The maximally extended spacetime continues indefinitely in all directions. 

The region covered by the cosmological coordinates (r, r) lies inside the bold lines. 
Two horizons are labeled. The inner horizon is the extension of the line labeled 
r = cxD to the region between the singularities (also see Fig. 5). The solid curves 
represent r = constant surfaces and the dotted curve shows a typical r = constant 
surface. The dashed fine denotes r — 0. 

Fig. 2: Penrose diagram for the Q — M RNdS geometry with M — Mext- The 
notation is the same as for Fig. 1. 

Fig. 3: Penrose diagram for the Q — M RNdS geometry with M > Mext- The 
notation is the same as for Fig. 1. 
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Fig. 4: Qualitative representation of the geometry on spacelike surfaces of 
constant r: (a) r < 0, (b) r = 0, (c) < r <^ Mi/ Ha, where a is the typical 
coordinate distance between the centers, (d) Mi/ Ha <^ r. 

Fig. 5: Regions of trapped surfaces and the horizons that separate them for 
the RNdS geometry. 

Fig. 6: A plot of the function F[R) = R + H{R + M)^. In determining the 
motion of radial null geodesies, one is lead to inequalities of the form dR/dy < F{R) 
or dR/dy > F{R). In the former case, it is clear that if it! < Ri,h initially, it will 
continue to decrease and become negative. These curves must hit the singularity. 
In the latter case, if it^^/j < R < R^s initially, then it must remain positive. These 
curves reach the de Sitter horizon. If equality holds, then R approaches Rds- 

Fig. 7: The motion of a Q = M shell in a RNdS geometry, (a) undermassive, 
(b) overmassive case. The shell follows a curve r = constant = tq. In the exterior 
region, r > ro (unshaded part of the diagram) the geometry is RNdS. In the 
interior (shaded, including the origin r — 0) spacetime is homogeneous de Sitter 
space. The shell eliminates the singularity that would otherwise be present in the 
right half of the digram (cf. Figs 1 and 3). Shells that replace the singularity on 
the left by a de Sitter interior follow a curve r' = constant. The corresponding 
Penrose diagrams would be reflections of (a) and (b) about a vertical axis. 

Fig. 8: (a) Penrose diagram for a two dimensional subspace of the two black 
hole solution. Each hole is undermassive, but their sum is overmassive. Near 
r — oo, the left part of the diagram is similar to Fig. 3 (single mass with M > Mext), 
and near r = the right half resembles Fig. 1 (single mass with M < Mext)- (b) 
Similar to (a) but with a shell replacing the singularity on the left. 
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